Parametric Oscillatory Instability in Fabry-Perot (FP) Interferometer 

V. B. Braginsky, S. E. Strigin and S. P. Vyatchanin 

Physics Faculty, Moscow State University, Moscow 119899, Russia 

e-mail: vyat@hbar.phys.msu.su 

(February 7, 2008) 



O 
O 

(N 

&q: 

< 

(N 
> 

On 

o 
r-» 

o 






% 



We present an approximate analysis of a nonlinear ef- 
fect of parametric oscillatory instability in FP interferom- 
eter. The basis for this effect is the excitation of the addi- 
tional (Stokes) optical mode with frequency Wi and of the 
mirror's elastic mode with frequency u m when the optical 
energy stored in the main FP resonator mode with fre- 
quency ojo exceeds the certain threshold and the frequen- 
cies are related as u>o — u>i +w m . This effect is undesirable 
in laser gravitational wave antennae because it may cre- 
ate a specific upper limit for the value of energy stored in 
FP resonator. In order to avoid it the detailed analysis of 
the mirror's elastic modes and FP resonator optical modes 
structure is necessary. 



I. INTRODUCTION 

The full scale terrestrial gravitational wave anten- 
nae are in process of assembling and tuning at present. 
One of these antennae (LIGO-I project) sensitivity ex- 
pressed in terms of the metric perturbation amplitude 
is projected to achieve soon the level of h ~ 1 x 10 -21 
PP]- I n 2008 the projected level of sensitivity has 
to be not less than h ~ 1 x 10~ 22 ||. This value 
is scheduled to achieve by substantial improvement 
of the test masses (mirrors in the big FP resonator) 
isolation from different sources of noises and by in- 
creasing the optical readout system sensitivity. This 
increase is expected to be obtained by rising the value 
of optical energy £q stored in the FP resonator opti- 
cal mode: £q > 30 J (it corresponds to the circulating 
power W bigger than 1 megaWatt). So high values 
of £q and W may be a source of the nonlinear ef- 
fects which will prevent from reaching the projected 
sensitivity of h ~ 1 x 10~ 22 . Authors of this article al- 
ready described two such effects: photo-thermal shot 
noise Q (the random absorption of optical photons in 
the surface layer of the mirror causes the fluctuating 
of mirror surface due to nonzero coefficient of ther- 
mal expansion) and photo-refractive shot noise H (the 
same random absorption of optical photons causes the 
fluctuations of the reflected wave phase due to the de- 
pendence of refraction index on temperature) . In this 
paper we analyze undesirable effect of parametric in- 
stability — another " trap" of pure dynamical nonlin- 
ear origin which (being ignored) may cause very sub- 
stantial decrease of the antennae sensitivity and even 
may make the antenna unable to work properly. 

It is appropriate to remind that nonlinear coupling 
of elastic and light waves in continuous media pro- 
duces Mandelstam-Brillouin scattering. It is a clas- 



sical parametric effect, however, it is often explained 
in terms of quantum physics: one quantum htoo of 
main optical wave transforms into two, i. e. %u>\ in 
the additional optical wave (Stokes wave: io\ < luq) 
and huj m in the elastic wave so that coq — U)i + <^m 
(it is Manley-Rowe condition for parametric process) . 
The irradiation into the anti-Stokes wave is also pos- 
sible (lui = ujq + uj m ), however, in this case the part 
of energy is taken from the elastic wave. The phys- 
ical "mechanism" of this coupling is the dependence 
of refractive index on density which is modulated by 
elastic waves. If the main wave power is large enough 
the stimulated scattering will take place, the ampli- 
tudes of elastic and Stokes waves will increase sub- 
stantially. The physical description is the following: 
the flux of energy into these waves is so large that be- 
fore being irradiated from the volume of interaction, 
the oscillations with frequencies io\ and to m stimulate 
each other substantially increasing the power taken 
from the main wave. Note that stimulated scattering 
causes irradiation only into Stokes wave because the 
additional energy pump into elastic wave must take 
place for radiation into anti-Stokes wave. 

In gravitational wave antennae elastic oscillations in 
FP resonator mirrors will interact with optical ones 
being coupled parametrically due to the boundary 
conditions on one hand, and due to the ponderomotive 
force on the other hand. Two optical modes may play 
roles of the main and Stokes waves. High quality fac- 
tors of these modes and of the elastic one will increase 
the effectiveness of the interaction between them and 
may give birth to the parametric oscillatory instability 
which is similar to stimulated Mandelstam-Brillouin 
effect H . This instability may create a specific upper 
limit for the value of energy £q . 

It is worth to note that this effect of parametric 
instability is a particular case of the more general 
phenomenon related to the dynamical back action of 
parametric displacement meter on mechanical oscil- 
lator or free mass. This dynamical back action was 
analyzed and observed more than 30 years ago (TJJ^]. 
Usually parametric meter consists of e.m. resonator 
with high quality factor Q (radiofrequency, microwave 
or optical ones) and high frequency stability pumping 
self-sustained oscillator. The displacement of the res- 
onator movable element modulates its eigenfrequency 
which in its turn produces the modulation of the e.m. 
oscillations amplitude (it can produce also phase or 
output power modulations). If the experimentalist 
attaches a probe mass to the movable element of the 
meter he is inevitably confronted with the effect of 
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dynamical back action, the ponderomotive force pro- 
duces a rigidity and due to finite e.m. relaxation time 
- a mechanical friction. Both these values may be 
positive and negative ones. In the case when the neg- 
ative friction is sufficiently high the behavior of me- 
chanical oscillator and meter becomes oscillatory un- 
stable. This effect was observed and explained for 
the case when value of mechanical frequency io m was 
substantially smaller than the bandwidth of e.m. res- 
onator MM ■ In this article we analyze the parametric 
oscillatory instability in two optical modes of FP res- 
onator and elastic mode of the mirror. In this case 
the value of u; m is much larger than the bandwidth of 
the optical modes. 

In section II we present the analysis of this ef- 
fect for simplified one-dimensional model which per- 
mits to obtain approximate estimates for the insta- 
bility conditions. In section III we present consid- 
erations and preliminary estimates for nonsimplified 
three-dimensional model, and in section IV — the pro- 



gram of necessary mode numerical analysis. 



II. SIMPLIFIED ONE-DIMENSIONAL MODEL 
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FIG. 1. Scheme of FP resonator with movable mirror 
(a) and frequency diagram (b). 

For approximate estimates we present in this section 
the simplified model analysis where we assume that: 

• The mechanical oscillator (model of mirror) is 
a lumped one with single mechanical degree of 
freedom (eigenfrequency io m and quality factor 

(°C m ^771/ ^0 771 ) ■ 

• This oscillator mass m is the FP resonator right 
mirror (see fig. H) having ideal reflectivity and 
the value of m is of the order of the total mirror's 
mass. 

• The left mirror (through which FP resonator is 
pumped) has an infinite mass, no optical losses 



and finite transmittance T = 2i\ L j '(XqQ op t) (Ao 
is the optical wavelength, Q op t is the quality fac- 
tor, L is the distance between the mirrors). 

• We take into account only the main mode with 
frequency loq and relaxation rate <5o = ^o/2Qo 
and Stokes mode with lu\ and 5\ = luq/2Qi cor- 
respondingly (Qo and Qi are the quality fac- 
tors), U!q — Ull ~ L0 m . 

• Laser is pumping only the main mode which 
stored energy So is assumed to be a constant 
one (approximation of constant field). 

It is possible to calculate at what level of energy £q 
the Stokes mode and mechanical oscillator becomes 
unstable. The origin of this instability can be de- 
scribed qualitatively in the following way: small me- 
chanical oscillations with the resonance frequency io m 
modulate the distance L that causes the exitation of 
optical fields with frequencies LOo±uj m . Therefore, the 
Stokes mode amplitude will rise linearly in time if time 
interval is shorter than <5f 1 . The presence of two opti- 
cal fields with frequencies loq and io\ will produce the 
component of ponderomotive force (which is propor- 
tional to square of sum field) on difference frequency 
io§—io\. Thus this force will increase the initially small 
amplitude of mechanical oscillations. In other words, 
we have to use two equations for Stokes mode and me- 
chanical oscillator and find the conditions when this 
"feedback" prevails the damping which exists due to 
the finite values of Q m and Q%. Below we present only 
the scheme of calculations (see details in Appendix |A|) . 
We write down the field components of optical 
modes and the displacement x of mechanical oscillator 
in rotating wave approximation: 

E Q = A [D e- iu,ot + D*e luJot ], 



* iu)±t] 



E x = ^^e'" 1 ' + D\e 
x = Xe~ iUmt + X*e iUmt , 

where Dq and D\ are the slowly changing complex 
amplitudes of the main and Stokes modes correspond- 
ingly and X is the slowly changing complex amplitude 
of mechanical displacement. Normalizing constants 
Ao, A\ are chosen so that energies So, l stored in each 



mode are equal to So l 



y 0,l 



Do,i| 2 /2. Then it is 



easy to obtain the equations for slowly changing am- 
plitudes: 



d t D 1 + 5 l D l = 



d t X + S m X = 



iX" 






i DqDIujqlji 



iAuit 



(1) 

(2) 



rmo m L 

where Alo = loo — W\ — Lo m is the possible detuning. 
Remind that we assume Do as a constant. 

One can find the solutions of (pi 0) in the following 
form D x {t) = D ie ( A - Aw / 2 )*, X*\t) = A* e ( A+A "/ 2 )* 
and write down the characteristic equation. The para- 
metric oscillatory instability will appear if real part of 
one of the characteristic equation roots is positive. 
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In LIGO design the values of Sq and <5i are of the or- 
der of the bandwidth the gravitational burst spectrum 
is expected to lie in, i.e. ~ 2ir x 100 s _1 . On the other 
hand many efforts were made to reduce the value of 
S m to the lowest possible level and thus to decrease the 
threshold of sensitivity caused by Brownian noise. In 
existing today fused silica mirrors Q m ~ 10 6 — 2 x 10 7 
and even for u> m = 10 7 s" 1 the value of S m < 10 s _1 . 
Thus we can assume that 6 m <C Si and obtain the 
instability condition in simple form: 



n 



K Q 



Aw 
WiW„ 



ImlP-uiL 5\8 n 



>1, 



So 



2£qQiQ 7 
mLP-iJL 



(3) 



(4) 



For estimates we assume parameters corresponding 
to LIGO-II to be: 



uj m = 2 x 10 5 sec 1 , 

<5i = 6 x 10 2 sec" 1 , 

£ - 3 x 10 8 erg, 

m — 10 4 g, 



sec 



L 



= 5 x 10" 



~ 2 x 10 if> sec 
= 4 x 10 5 cm 



(5) 



The mechanical frequency io m is about the fre- 
quency of the lowest mirror elastic (longitudinal or 
drum) mode and it has the same order as the inter- 
modal interval ~ irc/L ~ 2x 10 5 sec -1 between optical 
modes of FP resonator. The mechanical relaxation 
rate S m corresponds to the loss angle <f> ~ 5 x 10~ 8 
(quality factor Q m ~ 2 x 10 7 ) for fused silica. The 
value of energy £ corresponds to the value of cir- 
culating power about W ~ c£o/2L ~ 10 13 erg/s = 
10 6 Watt. 

For these parameters we have obtained the estimate 
of coefficient IZo for the resonance case (|Aw| <C <$i): 

TZo ~ 300 > 1 

It means that the critical value of stored energy £q 
for the instability initiation will be 300 times smaller 
than the planned value ~ 3 x 10 8 erg = 30 J q 

For nonresonance case and planned value of £ o the 
"borders" of detuning Auj cr it within the system is 
unstable, are relatively large: Ato cr it — Siy/Ro ~ 
1.7 x 10 4 sec" 1 . 



III. CONSIDERATIONS OF 
THREE-DIMENSIONAL MODES ANALYSIS 

The numerical estimates for the values of factor TZq 
and detuning Aw obtained in the preceding section 



It is worth to note that if sapphire is chosen then due to 
the larger value of Q m the factor 7?.o will be even bigger: 
7?.o — 5 x 10 3 . It is another argument which is not in favor 
of this material for mirrors. 



have to be regarded as some kind of warning about 
the reality of the undesirable parametric instability 
effect. In the simplified analysis we have ignored the 
nonuniform distribution of optical fields and of me- 
chanical displacements over the mirror's surface. It is 
evident that more accurate analysis has to be done. 
Below we present several considerations about further 
necessary analysis. 



A. The frequency range of "dangerous" optical 
and elastic modes 

The values of the mirror's radius R and thickness 
H for LIGO-II are not yet finally defined. Due to the 
necessity to decrease the level of thermoelastic and 
thermorefractive noises P,p|,p|pT|] the size of the light 
spot on the mirror's surface is likely to be substantially 
larger than in LIGO-I and the light density distribu- 
tion in the spot is not likely to be a gaussian one (to 
evade substantial diffractional losses) flO]j , Thus the 
presented below estimates for gaussian optical modes 
may be regarded only as the first approximation in 
which the use of analytical calculations is still possi- 
ble. 

The resonance conditions \ujo~ Wi— u m \ < Si maybe 
obtained with a relatively high probability for many 
optical Stokes and mirror elastic modes combinations. 
If we assume the main optical mode to be gaussian 
one with waist radius wq of the caustic (the optical 
field amplitude distribution in the middle between the 
mirrors is ~ e~ r ^ w °)), and if we assume also that 
the Stokes mode may be described by generalized La- 
guerre functions (Gauss-Laguerre beams) then the set 
of frequency distances Aw op t between the main and 
Stokes modes is determined by three integer numbers: 



7tc / 2(2N + M) LX Q 

Aw ptic ~ — I it H arctan ? 

L \ -k 2-kwq 



(G) 



where Ao is the wave length, K = Oil, ±2 ... is 
the longitudinal index, N = 0, 1, 2..., and M = 
0, 1, 2 . . . are the radial and angular indices. 

For u»o — 5.9 cm the beam radius on the mirror's 
surface is equal to w ~ 6 cm, corresponding to the 
level of diffractional losses about 20 ppm for mirror 
radius of R — 14 cm. In this case the equation (pi) has 
the following form: 



Auj. 



optic 



(2.4 if + 0.56 N + 0.28 M) x 10 5 s _1 . (7) 



We see that the distance between optical modes is not 
so large, i.e. ~ 3 x 10 4 s~ 4 . In units of optical modes 
bandwidth 2S X ~ 10 3 s^ 1 it is about 3 x 10 4 /2<5i ~ 30. 
Thus assuming that the value of elastic mode fre- 
quency can be an arbitrary one we can roughly es- 
timate the probability that the resonance condition is 
fulfilled as - 1/30. 

The order of the distance Aw m o between the fre- 
quencies for the first several elastic modes is about 
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Aoj m o ~ Trv s /d ~ 2 x 10 5 s" 1 (d is the dimension of 
the mirror and v s is the sound velocity). It is about 
one order larger than the distance between the opti- 
cal modes. However, for higher frequencies uj m these 
intervals become smaller and can be estimated by for- 
mula 



Aw„ 



2ojI 



Even for io m ~ 6 x 10 5 s 1 the intervals between the 
elastic and optical modes become equal to each other 



and has value about ~ 3 x 10 s 



And for 



10 7 



the distances between elastic modes become 



of the order of optical bandwidth 2<5i . Therefore, the 
resonance condition for these frequencies is practically 
always fulfilled. 

On the other hand according to (Q) the factor 7?.o 
decreases for higher elastic frequencies uj m . In addi- 
tion the loss angle in fused silica usually slightly in- 
creases for higher frequencies [0,|l2|. Assuming that 
the upper value cu m ~ 2 x 10 6 s" 1 , Q m ~ 3 x 10 6 and 
other parameters correspond to (g) we obtain IZq ~ 1 . 
Therefore, the elastic modes which " deserve" accurate 
calculations lie within the range between several tens 
and several hundreds kiloHerz. The total number of 
these modes is about several hundreds. 



B. The Matching between the Mechanical 
Displacements and Light Density Distributions 

The simplified model described in section O is ap- 
proximately valid for the uniform over all the mir- 
ror's surface distribution of optical field density and 
pure longitudinal elastic mode. The equations for this 
model can be extended for any distributions of me- 
chanical displacements in the chosen elastic mode and 
for any distribution of the light field density in chosen 
optical modes. This extension can be done by adding 
a dimensionless factor A in (ph: 



^oA 



Su 2 



>1, 



A 



^(//oCrxJ/iCfxKrfn.) 2 
t\h?dr±j\h\ 2 dr ± j\v\HV 



(8) 
(9) 



Here /o and f\ are the functions of the distributions 
over the mirror's surface of the optical fields in the 
main and Stokes optical modes correspondingly, vec- 
tor u is the spatial vector of displacements in elastic 
mode, u z is the component of u, normal to the mirror's 
surface, J dr± corresponds to the integration over the 
mirror's surface and J dV — over the mirror's volume 
V. 

It is necessary to know the functions /o , /i , u in or- 
der to calculate the factors A for different mode com- 
binations. But to the best of our knowledge there is 
no analytical form for u in the case of cylinder with 



free boundary conditions [JO]]. Our approximate esti- 
mates show that there is substantial number of modes 
combinations for which factor A is large enough to 
satisfy the condition (g) for parameters (|5|). We do 
not present the details of these estimates here because 
they are rather rough. It is evident that a complete 
numerical analysis which includes nongaussian distri- 
butions of optical fields is necessary to do. 



IV. CONCLUSION 

The simplified model analysis of parametric oscil- 
latory instability and considerations about the real 
model presented above may be regarded only as the 
first step along the route to obtain a guarantee to 
evade this nondesirable effect. Summing up we may 
formulate several recommendations for the next steps: 

1 . Due to the finite size of the mirror and to the use 
of the nongaussian distribution of light density 
we think that the accurate numerical analysis 
of different optical and elastic mode combina- 
tions (candidates for the parametric instability) 
is inevitably necessary. This problem (numeri- 
cal calculations for the elastic modes) has been 
already solved partialy pi . 

2. In the same time the numerical analysis may not 
give an absolute guarantee because the fused sil- 
ica pins and fibers will be attached to the mirror. 
This attachment will change the elastic modes 
frequency values (and may be also the distribu- 
tion) . In addition the unknown Young modulus 
and fused silica density inhomogeneity will limit 
the numerical analysis accuracy. Thus the di- 
rect measurements for several hundreds of probe 
mass elastic modes eigenfrequencies values and 
quality factors are also necessary. 

3. When more "dangerous" candidates of elastic 
and Stokes modes will be known, their undesir- 
able influence can be depressed. For example 
it can be done by small the change of mirror's 
shape. 

4. It is also reasonable to perform direct tests of 
the optical field behavior with smooth increase 
of the input optical power: it will be possible 
to register the appearence of the photons at the 
Stokes modes and the rise of the Q m in the corre- 
sponding elastic mode until the power W in the 
main optical mode is below the critical value. 

5. Apart from above presented case of the oscilla- 
tory instability it is likely that there are similar 
instability in which other mechanical modes are 
involved (especialy violin ones which also have 
eigen frequencies several tens kHz and higher). 
There are also additional instability for the pen- 
dclum mode in the mirror's suspension (in the 
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case of small detuning of pumping optical fre- 
quency out of resonance) . These potencial "dan- 
gers" also deserve accurate analysis. 

We think that the parametric oscillatory instabil- 
ity effect can be excluded in the laser gravitational 
antennae after this detailed investigation. 
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APPENDIX A: LAGRANGIAN APPROACH 

Let us denote qq (t) and q\ (t) as generalized coordi- 
nates for the FP resonator optical modes with frequen- 
cies wo and u>i correspondingly, so that their vector 
potentials (A , A^, electrical (E , E{) and magnetic 
(Hq, Hi) fields are the following: 



Mt) 



I2nc 2 



(ti 



tie 



-ik-i 



Qi(t), 



2tt 



Ei(t) = -\jj£ (he lk > z - f*e- ikiZ ) dtq.it), 



H t (t) 



2tt 



Q r \fi 



ik t z 



* —ikiz\ 



SiL 

fi = fi(r±,z), Si 



tie 



tOiqi{t), 



\M 2 dr ± . 



Let also denote x(t) as generalized coordinate of the 
considered elastic oscillations mode with displacement 
spatial distribution described by the vector u(r ) . Now 
we can write down the lagrangian: 



C — Cq + JL-x + C r , 



*-*?' 



/;i = L((E Q ) 2 - (H ) 2 ) & ^ _d t q 2 o;ggg 



Ci = 

r — 



8tt 

dtgj _ iolgj 
2 2' 

M{d t x) 2 Muj 2 m x 2 



We consider only one mechanical mode below. Now 
we can write down the equations of motion (adding 
losses in each degree of freedom): 

d 2 qo + 25od t qo + ^qo = -By WoWigi, 

2x 
d 2 q x + 25id t qi + to\qi = -B — ujiWo go, 



d 2 x + 26 m d t x + u 2 n x = -B 



L 

2cl>oCl>i 

ML 



qoQi- 



Introducing slowly varying amplitudes we can rewrite 
these equations as: 



d t Do 



q {t)=D {t)e- iu ° t 

q 1 (t)=D 1 (t)e-^ t 
x(t) =X(t)e- luJmt - 

Aw = LUo — Wl — L0 m , 

iBXDiUi iAujt 
Oo-^o — ; e , 



-D*(t)e luJot , 

-Dl{t)e Wl \ 
X*(t)e iUmt , 



d t D 1 + 6 1 D 1 = iBX * DoUJo e-^ t 



d t X 



L 
5 m X = tggggigj* e -iA. t> 



(Al) 
(A2) 



We can see that this system (A2, Al) coincides with 
(|0) if B = 1, and A/ = to. 

For the simplest re sona nce case 8ui — it is easy to 
substitute (A2) into (|Al[ ) and to obtain the condition 
of parametric instability (in the frequency domain) : 



Di(b\ — til) = x 



-i BDqDiuiquji 



1 < 



L Mw m L(5 m + iQ) ' 

Condition of instability: 

£ 2 |A)| 2 w V 

Muj m L 2 Si5 m 



Now we can express the energy £q in mode "0" in 
terms of |-Do| 2; 



£o = 



dtql 



2 2 

^0% 



= - (i-iLOof [Do e-^ ot - D* e^ ot ] + 

+ io 2 [D e-^ 0t +D*e luJ0t ] 2 ^ = 

= 2lu 2 \D q \ 2 . 

Now we can write down the condition of parametric 
instability in the following form: 



M = p \u(r)\ 2 dV, 
Jv 

xu^Ho + Hrf 



*->int — 



dr i = 



B 



= -2uJooJ\qoq\B —, 

J fo(r±)fi(r ± )u z drj_ 

^I\M 2 df ± J\h\ 2 dn 



SoB 



W]W„ 



>1, 



A 



2m<J 2 n L 2 SiS n 

B 2 m = _ V (J fo(r±)fi{r±)u z df A 
M 



! \h\ 2 d^ $ \h\ 2 dr^ $ \u\ 2 dV 



(A3) 
(A4) 



which accurately coincides with (§L g). 

Let us deduce the instability condition for nonreso- 
nance case. We are looking for the solution of (EI, 0) 
in the following form: 
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G 



A_ = A 



D 1 {t)=D l e 

iAuj 



A_t 



A, 



A 



X*(t) = X*e A +\ 
iAui 



2 ' + 2 ' 

and writing down the characteristic equation as: 

DqWqUiA 



(X+ + S 1 )(X_ + 6 m )-A = 0, 
The solutions of characteristic equation are: 

61 + s„ 



A. 



Al,2 — 

Dot = 



2 
Si - S„ 



± VDet, 



iAu 



A. 



2 2 

The condition of instability is the following: 

£1 + S m 



WDet > 



Using a convenient formula: 



(A5) 



SRVDct 



Dct = a + ib, 
V2 



\fVa 2 + b 2 + ft. 



we can rewrite the condition (A5) as: 

2 



a 2 + b 2 + a) > 



Si + 5„ 



A 2 



-2A 



(Sj-Sm) 2 Au 2 

4 4 

(Si - 8 m ) 2 Alu 2 



(A6) 



(A7) 



No te th at for t he resonance case (Auj = 0) the solution 
of (A5 or A6) is known: A > SiS m . For our case 
S m "^ <5i it means that A <C S 2 . Therefore for small 
detuning Au> <C S± we can expand a 2 + b 2 in series in 
terms of A and rewrite condition (A6) as: 



!(v?TF + .K + & 



4 

A (SlS m ) 2 _ Auj 2 
2_ 4 4_ 

4 ^4 



A 4 



(&i-S m ) 2 , Ac^ 

4 '4 



> S^r, 



Or A > (5i<5 m x 



(Si 



Ao 



(*i - s m y 



(A8) 



Let us underline that condition (A8) is obtained for 
small detuning Aw <^ 6\. However, considering situa- 
tion more attentively one can conclude t hat expansion 
in series (and consequently the formula (A8) ) is valid 
for the condition: 



A« 



(6 1 -6 m ) 2 , Aco 2 



(A9) 



We see that this condition is fulfilled for the solution 
(A8). Therefore we conclude that solution (^ 



is 



approximately valid for any detunings Auj. 
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